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A COMPUTER SIMULATION OF THE AFTERBURNING PROCESSES OCCURRING WITHIN 


SOLID ROCKET MOTOR PLUMES IN THE TROPOSPHERE 

Richard I. Gomberg and Roger B. Stewart 
Langley Research Center 


SUMMARY 

As part of a continuing study of the environmental effects of solid rocket 
motor (SRM) operations in the troposphere, a numerical model is used to simulate 
the afterburning processes occurring in solid rocket motor plumes and to predict 
the quantities of potentially harmful chemical species which are created. The 
calculations include the effects of finite-rate chemistry and turbulent mixing. 
It is found that the amount of NO produced is much less than the amount of HC1 
present in the plume, that chlorine will appear predominantly in the form of HC1 
although some molecular chlorine is present, and that combustion is complete as 
is evident from the predominance of carbon dioxide over carbon monoxide. 


INTRODUCTION 

This is the second in a series of papers dealing with chemical species 
released into the troposphere as a result of solid rocket motor operations. 
Reference 1 (NASA TN D-8137) dealt primarily with NO production from the space 
shuttle and Titan III-C solid motors. The present study expands the scope of 
investigation to include the partitioning of chlorine between the compounds HC1, 
Cl, Cl 2 , and CIO. In addition, the partitioning of carbon between CO and C0 2 is 
considered, as well as the production of NO and N0 2 . 

The primary analytical tool used in the investigation was the low-altitude 
plume program, the LAPP code. (See ref. 2.) This computer program models the 
plume as axisymmetric , reacting concentric streamtubes with turbulent mixing 
between tubes. It is suitable for relatively low-altitude computational stud- 
ies of afterburning in rocket plumes. Above an altitude of 15 km, the rocket 
nozzle flow becomes increasingly underexpanded and a shock pattern develops in 
the plume that is not accounted for by the code. The section that follows 
describes how the LAPP program is used in the present study as well as a modifi- 
cation which can be used to strengthen weaknesses in the model. Subsequent sec- 
tions give estimates of the nitrogen, chlorine, and carbon-containing species 
which result from afterburning at various altitudes. 


SYMBOLS 


A 


area of plume, m 2 


B 


activation energy 


M Mach number 

n real number 

p pressure 

R gas constant 

r effective nozzle-exit radius (derived from true exit radius by means 

of an isentropic expansion) 

T temperature 

U velocity of particle within plume 

U av average velocity of particles within plume approaches 

x distance downstream of nozzle-exit plane, m 

y distance measured from center line of plume perpendicular to down- 

stream direction 

p density of plume 

Subscripts : 

e pertaining to exit plane of nozzle 

A^O^ pertaining to 

x pertaining to a distance x (meters) downstream of exit plane 

ANALYTICAL TECHNIQUES 
The LAPP Program 

Figure 1 is a diagram of a typical low-altitude rocket plume. The LAPP 
code used to model this plume is a time-independent axisymmetric model that 
divides the hydrodynamic flow field into streamtubes. The model is formulated 
in streamtube coordinates. The unperturbed atmosphere or free-stream flow is 
the outermost streamtube; the remaining streamtubes contain the plume. In each 
streamtube the governing differential equations are the reacting boundary-layer 
equations. (See ref. 3 .) The numerical scheme used to solve these equations 
is mixed explicit-implicit. 

Several eddy-viscosity models are presently available which describe the 
turbulent entrainment of air into the plume. The Donaldson-Gray model (ref. 4) 
was used in this investigation because of its success in prior studies by other 
investigators (ref. 5). 
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As inputs to the code, it is necessary to specify the conditions at the 
nozzle-exit plane (such as velocity, temperature, and composition of the exhaust 
gases), conditions in the unperturbed atmosphere, the chemical reactions to be 
considered, and all the pertinent thermochemical data. The nozzle-exit plane 
conditions for the space shuttle solid rocket motor (SRM) were furnished by 
Marshall Space Flight Center and the data for the Titan III-C SRM by Aerospace 
Corporation. (See tables I and II.) The unperturbed atmospheric conditions 
were taken from U.S. Standard Atmosphere (ref. 6); thermochemical data were 
taken from the JANNAF reports (ref. 7); and chemical reaction rate coefficients 
were taken from references 8 to 13. 

Calculations were made for the shuttle SRM at altitudes of 0.7, 6, 10, and 
15 km and for the Titan at altitudes of 0.9 and 18 km. At high altitudes the 
pressure of the gases at the exit plane is on the order of 1 atmosphere , whereas 
the ambient pressure is about 0.1 atmosphere. (1 atmosphere = 101.3 kPa.) For 
these high altitudes, a hand-calculated isentropic expansion is used to gain 
effective nozzle-exit plane conditions. This approach was used for the shuttle 
motor at 6, 10, and 15 km and for the Titan at 18 km. The chemical composition 
at the effective nozzle-exit plane was assumed to be identical to those of 
tables I and II. By using the mass flow rates taken from reference 14 and the 
isentropic assumption, the complete set of effective nozzle-exit plane condi- 
tions necessary to start the code at these higher altitudes are given in 
table III. 

Because the reaction mechanism used in the code is not appropriate for tem- 
peratures below 600 K, the code cuts off all chemical kinetic reactions below 
this temperature. At all altitudes between sea level and 15 km, both the shuttle 
and Titan motors produce regions extending at least 0.5 km downstream of the 
nozzle-exit plane in which temperatures are well above 600 K, 

The chemical reactions and the rate coefficients used in this study are 
given in table IV. These reactions are descriptive of processes occurring above 
600 K. Although they are chemically inert in this model, certain species such 
as Al^Oo are considered because of their importance in maintaining a heat bal- 
ance. These species are listed at the end of table IV. All constituents of the 
plume are taken to be ideal gases . 


Mass Conservation 


As was shown in reference 1, the LAPP code does not conserve mass. This 
effect is demonstrated in figure 2, which plots the mass flow of the inert spe- 
cies Al^Oo against downstream distance in the plume of a shuttle solid rocket 
motor at an altitude of 15 km. Since Al o 0 o is inert and the code is time- 


independent, this quantity should be conserved. If 
stream, the mass flow rate of A1 2 0^ is (pUA) A1 Q , 


x is the distance down- 
This quantity should be 


constant for all x. As can be seen from figure 2, the quantity actually 
decreases with downstream distance . 


Several calculations with variation in step size demonstrated that the 
problem is not a straightforward roundoff error as the mass loss was insensi- 
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tive to the step size used . The problem appears to be connected with the way 
the model allows the plume to spread out . In LAPP , the outermost streamtube is 
the unperturbed atmosphere which is taken to be a mixture of N2 and CU only. At 
the atmospheric boundary the code truncates any constituents in the plume such 
as AlgO,, HC1 , NO, etc. Thus, if the computational mesh is not allowed to 
spread as rapidly as the mass is viscously mixing outward, mass will in effect 
be lost by leaking out of the edges of the flow field . 

The mechanism for expanding the plume in the LAPP code is the addition of 
streamtubes. Thus, it is important to add streamtubes at a proper rate to 
account for the expansion of the plume and to satisfy the nonleakage boundary 
condition . In the unmodified version of the code , streamtubes are added when- 
ever the outermost tube (the unperturbed atmosphere) and its adjacent tube sat- 
isfy either 


u outer " U adjacent > Q 
^outer I 


or 


n - t 
outer adjacent 

T 

A outer 


> 0.05 


By varying the numbers which the inequalities must satisfy, one can control 
the rate of expansion of the flow field. At distances close to the nozzle-exit 
plane , the inequalities as they are stated seem to be adequate to satisfy the 
necessary boundary conditions. At distances far downstream, roughly greater than 
250 m, the inequalities are too steep to allow the plume to expand properly in 
this region . A modification that can be made , then , is to lower the inequalities 
necessary to add streamtubes far downstream. The amount of lowering is easily 
determined by keeping careful track of pUA | UA * Al 2 0g j • By 

choosing inequalities so that the resultant expansion causes this ratio to stay 
close to one, the model will satisfy the nonleakage boundary condition. 


Mass Conservation in Present Studies of Solid Booster Motor Plumes 

By making use of the modified criteria stated previously, the mass loss due 
to "leakage" can be held to a few percent. Consider, for example, figure 2(a) 
in which is plotted the mass flow of Al^O^ against downstream distance for a 
shuttle SRM at an altitude of 15 km. As can be seen in contrast to that behav- 
ior in figure 2(b) by carefully controlling the expansion of the plume, the mass 
flow oscillates with small amplitude about the conserved value. In figure 3, 
the mass flow of A^O, for a shuttle motor at an altitude of 5 km is plotted as 
a function of downstream distance. Again, the mass loss is seen to be small. 

For each case studied, a careful analysis of the mass flow rate of inert species 
using the techniques described here was carried out in order to investigate pos- 
sible mass losses. All calculations performed in this study were taken to a 
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distance of 1 km downstream of the exit plane . A reliable indicator of how 
closely mass is being conserved is obtained by comparing the mass flow of Al 2 0g 
at a distance of 1 km downstream of the exit plane with the mass flow at the 
exit plane. In table V such a comparison is made for the various cases already 
described . 


A Sensitivity Study 

In addition to calculating the behavior of the afterburning plume for the 
cases already described, an attempt was made to determine the sensitivity of the 
model to variations in the input parameters. In particular, chemical reaction 
rates are not known precisely, but have associated with them uncertainties which 
can be as high as a factor of a thousand. Eddy diffusion coefficients are also 
quantities which are not known with great precision. A variation in a rate 
coefficient will cause variation of the resultant chemical composition of the 
plume. However, if the change in composition is small even when the rate coef- 
ficient is varied over its full range of experimental' uncertainties, then the 
composition is said to be insensitive to uncertainties in the rate coefficient. 

To understand how sensitive this particular calculation is to experimental uncer- 
tainties in the rate coefficients, a series of computer simulations with varied 
rate coefficients were run of the afterburning from a shuttle SRM at an altitude 
of 10 km. 

In each of the cases in this sensitivity study, identical nozzle-exit plane 
conditions appropriate to the shuttle SRM at 10 km were used. The chemical 
scheme, however, is changed several times to investigate the uncertainties in 
rate coefficients in some important reactions. Consider, for example, table VI. 
Of the 36 reactions used in the chemical scheme, 12 directly include the chlo- 
rine atom. The forward kinetic rates associated with these reactions are made 
to vary in one simulation as shown in the first column of table VI. All rate 
coefficients used in this model are of the form 


Rate coefficient = CT n exp — 

RT 

where C is a constant and B is the activation energy. In the first column 
of table VI t new values for C are used for the 12 reactions shown. In the 
simulation, only these rate coefficients are varied; all other rates remain as 
shown in table IV . 

In the next simulation, the rate coefficients shown in the second column 
were used. These variations cause all reactions producing atomic chlorine to be 
accelerated (first column) in the first simulation and slowed (second column) in 
the next simulation. 

In table VII a similar study is made for atomic oxygen. One set of rate 
coefficients enhances 0 production, another set suppresses 0 production. In 
table VIII a set of rate coefficients speeding the combustion processes (within 
experimental uncertainties, of course) are used. In all cases, only the rate 
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coefficients shown are varied. As will be shown in detail, the model is fairly 
insensitive to the uncertainties in the chemical rate coefficients. 

To study the effect of variations in turbulent mixing, a case is run with 
a different mixing model . A good alternate to Donaldson-Gray is the Ting-Libby 
model (ref. 15) which is used in this last case of the sensitivity study. 

In each of the sensitivity study computer simulations, the conservation of 
mass was closely watched. As shown in table V, mass is conserved fairly accu- 
rately in these runs also . 


RESULTS AND DISCUSSION 

Calculations were made for a Titan solid propellant booster at 0.9 km and 
18 km, and a shuttle motor at 0.7, 6, 10, and 15 km with sensitivity calcula- 
tions at 10 km. These studies neglect the shock structure present in the plume 
and introduce some uncertainties, particularly at high altitudes. With increas- 
ing x/r, the temperatures in the plume decrease, the chemical reactions begin 
to slow down, and the velocities decrease. As x/r approaches the value where 
the center-line temperature drops below 600 K, chemical reactions are shut off 
in the model . Although the plume will continue to diffuse and therefore the 
absolute concentration of each specie will decrease with increasing x/r, the 
total mass flow of each specie crossing a plane perpendicular to x will be 
constant in this region (within the mass conservation limits previously 
described ) . 

Farther downstream, the average velocity of the particles in the plume 
approaches ambient velocity. (All calculations are carried out in the frame of 
reference of the nozzle exit plane.) Thus, not only is the mass flow (PU ay A) x 
constant, but also 


(p0 av A >x 

<“av>x 


(PA) x ► Constant 


The physical meaning of (pA) is that quantity of mass in a slab, the 
dimensions of which are equal to the area of the plume whose plane is perpendic- 
ular to x and whose length is 1 m in the x-direction. The value of (pA) x 
for a given species for large x is a reasonable indication of the quantity of 
that species that is deposited in the atmosphere. The unit of (pA) x is g/m. 


NO Production in Solid Booster Motor Plumes 

Table IX lists the production levels of NO and NO2 for the various cases 
previously described. There are a number of interesting features in these 
results. At altitudes below 6 km, the shuttle booster produces roughly twice 
the NO of the Titan booster. This is quite reasonable since the two motors are 
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similar and the mass flow of the shuttle is roughly twice that of the Titan. At 
the high altitude of 18 km, however, in spite of the larger mass flow of the 
shuttle, the Titan SRM produces roughly twice the NO of the shuttle at an alti- 
tude of 15 km. As discussed in reference 1 , this difference is due to the exis- 
tence of a thrust-vector control system present in the Titan, but not in the 
shuttle. The thrust-vector control system consists of a tank of N 2 0ij which is 
strapped onto each booster and which is bled into the chamber through a series 
of steering jets in order to make fine course adjustments. In the high tempera- 
tures present in the chamber , much of the N^O^ breaks down into NO as well as 0 
and 0 2 . ^ 

The effects of uncertainties in the rate coefficients on the production of 
NO are inherent in the results of the parametric studies of the shuttle at 10 km. 
The rate coefficients that are being used as best estimates result in a (PA) X 
of NO of 25 g/m at large x. However, when some of the uncertainties of the 
rate coefficients are put into the calculations, the result is that this NO pro- 
duction can be as low as a (PA) X 21 g/m or as high as 37 g/m for large x. 

By using the Ting-Libby viscosity model instead of the Donaldson-Gray model, 
the code predicts that oxygen is entrained into the plume at a more rapid rate 
than Donaldson-Gray predicts. This condition causes the combination of nitrogen 
and oxygen to occur at higher temperatures and over a longer period of time. Cor- 
respondingly, the NO levels rise to 41 g/m from 25 g/m. 

The NO levels can be seen in some detail in figure 4. The quantity of NO 
and various other species in the plume in units of g/m is plotted as a function 
of downstream distance for the shuttle motor at 10 km. Results for both the 
Donaldson-Gray and Ting-Libby models are shown. For both models the quantity of 
A1 2 0^ as a function of downstream distance is very close. For the Donaldson- 
Gray model the NO production can be seen to rise to an asymptotic value of 31 g/m 
whereas it rises to a value of 41 g/m for the Ting-Libby case. Note that for the 
Ting-Libby case, the NO levels rise more rapidly than for the Donaldson-Gray case 
in the early parts of the plume where T is higher. 


Chlorine Partitioning in Solid Booster Motor Plumes 

Table X shows how chlorine is partitioned among HC1 , Clg, and CIO in the 
various situations discussed. In all cases the plume partitions chlorine in 
amounts by weight of roughly 89 percent HC1, 9 percent Cl 2 , and about 2 percent 
is in the form of Cl or CIO. By mole fraction HC1 makes up more than 9 percent 
of the reacting chlorine-containing compounds in the plume. The shuttle plume 
will also contain some nonreacting particles such as A1C1 and A1C1 2 , 

This partitioning of the reacting chlorine compounds is true for both the 
Titan and shuttle motors. In figure 5, (PA) X of various chlorine-containing 
species are plotted against downstream distance for a Titan motor at 0.9 km. In 
figure 6 a similar plot is made for a shuttle motor at an altitude of 5 km. In 
both cases as is the case at other altitudes, the partitioning of chlorine fol- 
lows the pattern described . 
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Completeness of Combustion in Solid Booster Motor Plumes 


The completeness of combustion can be discerned from the ratio of C0 2 to 
CO. At the exit plane of the Titan, there is a 9 to 1 ratio in g/m of CO to g/m 
of C0 2 . For the shuttle motor there is a 7 to 1 ratio of CO to C0 2 at the exit 
plane. By the time afterburning is complete (at approximately 1 km downstream), 
much of the CO has combined with oxygen to form C0 2 . The extent to which C0 2 
has been created at the expense of CO is a good indication as to how complete 
the combustion is in the plume. From table XI it is clear that the combustion 
is almost total. Ratios of C0 2 to CO of a 100 to 1 are typical of the after- 
burned products. 

The process of converting CO to C0 2 for a Titan motor is shown in figures 7 
and 8. Figure 7 shows (PA) X in g/m of various species against downstream dis- 
tance in the plume of a Titan motor at 0.9 km. As can be seen from figure 7, at 
the exit plane CO is the predominant species, but as a result of afterburning, 
C0 2 becomes by far the more prevalent quantity by the time the asymptotic region 
has been reached. In figure 8 is plotted the mole fraction of several species 
along the center line of the plume against downstream distance, again for the 
Titan at 0.9 km. Here again, the conversion of CO to C0 2 can be seen as the CO 
cui?ve falls off rapidly when afterburning sets in and the mole fraction of C0 2 
increases at the expense of the amount of CO present . 

In figures 9 and 10 a similar set of curves is shown, for the shuttle at 
6 km. Both for (PA) X of species (fig. 9) and for the center-line mole frac- 
tion curves (fig. 10), the conversion of CO- to C0 2 is easily seen, 


Composition and Physical Characteristics 

Figures 11 to 20 present a detailed analysis of the composition and physi- 
cal characteristics of the plume of a shuttle motor at 0.7 km and at 15 km. Fig- 
ures 11 to 13 describe (PA) (in g/m). Notice that the (PA) X of "combustion" 
species produced at 0.7 km (fig. 11) are much greater than those at 15 km. This 
difference is due to the following condition: (1) the mass flow at the exit 

plane is smaller at 15 km than that at 0.7 km (ref. 14), (2) the free-stream 
velocity at 15 km is .much greater than that at 0.7 km, and (3) the free-stream 
pressure at 15 km is less than that at 0.7 km. 

In figure 12, quantities of chlorine-containing species (pA) are plot- 
ted as a function of downstream distance for the 0.7- and 1.5-km shuttle motors. 
These figures show the partitioning of chlorine among the various species, 
whereas figure 13 concentrates on oxides of nitrogen. 

Figures 14 to 16 show for the same conditions and species as figures 11 
to 13, the behavior of the mole fraction along the center line plotted against 
downstream distance . In these figures it is seen that the mole fractions along 
the center line of the 15-km plume follow a pattern similar to those at 0.7 km. 
However, at this high altitude, since there are fewer collisions/meter the pro- 
cesses are stretched out over greater downstream distances. 
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In figures 17 and 18 are cross-sectional velocity profiles for the 0.7- and 
15-km shuttle motors at several downstream stations. In figures 19 and 20 tem- 
perature profiles for these two cases are presented. At x/r = 0, both tempera- 
ture and velocity across the exit plane are assumed to be constant and are much 
higher than ambient levels. As the plume gases entrain air and afterburn farther 
downstream, the temperature particularly along the center line rises even above 
the nozzle-exit plane value. Upon completion of afterburning, as dilution sets 
in, the plume cools so that far downstream the temperatures of the plume gases 
approach ambient values. The velocity of the gases is a decreasing function of 
downstream distance as the shear forces slow the gas molecules to ambient speeds. 

Table XII contains a summary of the chemical composition of all plumes con- 
sidered in this study. 


CONCLUDING REMARKS 

The physical, thermodynamic, and chemical properties of solid rocket motor 
plumes have been studied by use of an axisymmetric , time-independent model which 
includes the effects of turbulent mixing and finite rate chemistry. A modifica- 
tion has been introduced into the model to conserve mass far downstream of the 
exit plane. 

The shuttle and Titan III-C boosters are studied. The quantity of NO pro- 
duced seems to fall off rapidly with increasing altitude. Combustion is seen to 
be complete as nearly all the CO at the exit plane has been converted to C0£ far 
downstream. Chlorine is partitioned largely between HC1 and Cl 2 with 90 percent 
of the weight in the form of HC1 . . 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 6, 1976 
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TABLE I.- SPACE SHUTTLE SOLID ROCKET MOTOR EXIT PLANE CONDITIONS 

[This table is based on an equilibrium model. It is appropriate 
for a shuttle motor at 0.7 km altitude. Space shuttle data 
were obtained from Ben Shackleford, Marshall Space Flight 
Center, Mar. 1974. These data are appropriate for altitudes 
less than 1 km. The mass flow obtained from using these 
inputs is 25 percent lower than that predicted in ref. 14.] 


Species 

Mole fractions 

A1C1 

0.00003 

aici 2 

.00008 

Aid, 

.00001 

A10CI 

•00002 

Al-Oo 

.07975 

CO 

.23205 

C0 o 

.02120 

Cl 2 

.00225 

Fe 

.00011 

FeCl 

•00001 

FeCl 0 

.00122 

H 

.00556 

HC 1 

.15442 

H 2 

.27812 

h 2 o 

. 14068 

NO 

.00001 

n 2 

.08401 

0 

.00001 

OH 

.00046 

CIO 

8.1 x nr 9 

Cl 2 

6.9 x 10 -7 

HCU 

1 .2 x 10 -10 

2 

no 2 

5.8 x lO " 1 ^ 

n 0 6 

8.6 x 10" 11 

^2 

Op 

2.0 x 10 -7 

N 2 

.0000 


Exit plane conditions 


p, atm . 

T, K . . 

M . . . 
r , m . . 

U, m/sec 


0.8345 
2308 
2.779 
1 .810 
2433.8 
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TABLE II.- TITAN MOTOR DATA 


[This table is based on an equilibrium model. It is appropriate for 
altitudes of approximately 1 km. In the third column, it is 
assumed that the injected NgOjj breaks down into NO and 02*] 


Exit plane 
species 


CO 
OH 
CO 2 
H 

NO 

no 2 

N 

Ih 

NpO 

Cl 
Cl 2 
CIO 
A1 2 0 3 


Without thrust vector 
control system 


2.468 

X 

10" 1 

2.452 

X 

10 

9.000 

X 

10~5 

9.000 

X 

10' 

1 .811 

X 

10“ 2 

1 .811 

X 

10 

3.200 

X 

10‘ 3 

3.200 

X 

10 

1 .035 

X 

10 -3 

1 .034 

X 

10 

5.858 

X 

10 _1 1 

5.585 

X 

10 

3.157 

X 

10" 1 

3-136 

X 

10 

1 .058 

X 

10“! 

1 .051 

X 

10 

9.000 

X 

10-5 

6.680 

X 

10 

2.002 

X 

10-1 1 

2.002 

X 

10 

1 .884 

X 

10 -8 

1 .884 

X 

10 

8.105 

X 

io~ 2 

8.051 

X 

10 

1 .510 

X 

10 ft 

1 .500 

X 

10 

1 .090 

X 

10“® 

1.090 

X 

10 

1 .150 

X 

10“ 3 

1 .150 

X 

10 

3-293 

X 

10 ‘ft 

3.293 

X 

10 

2.392 

X 

10 -8 

2.392 

X 

10 

7.680 

X 

W ~7 

7.170 

X 

10 

6.959 

X 

10" 7 

6.959 

X 

10 


With thrust vector 
control system 
(based on 
50 lb/sec-motor) 


)-•> 

r 5 

)- 2 

)-3 

)-3 

r 11 

I 

% 

$ 

>- 2 

i-7 


Exit plane conditions 

Mass flow, Mg/sec 1.84 

T e , K 1921 

P e * atm 0.702 

r e » m 1-35 
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TABLE III.- INITIAL CONDITIONS 


[By using these inputs, the predicted mass flows agree within 
10 percent with those in ref. 14 for the shuttle] 


(a) Ambient conditions 


Altitude , 
km 

Vehicle 

Ambient 
pressure , 
atm 

Ambient 
temperature , 
K 

Free-stream 

velocity, 

m/sec 

0.7 

Titan 

0.702 

300 

30.48 

6 

Titan 

• 536 

260 

275.8 

18 

Titan 

• 0797 

207 

537-6 

-9 

Shuttle 

.8345 

300 

30.48 

5 

Shuttle 

.536 

260 

259-3 

15 

Shuttle 

.119 

216 

552.2 


(b) Exit plane conditions 


Altitude , 
km 

Vehicle 

Effective exit 
radius , 
m 

Effective exit 
temperature , 

K 

Effective exit 
velocity, 
m/sec 

0.7 

Titan 

1 .35 

1921 

2501 

6 

Titan 

1 -35 

1921 

2501 

18 

Titan 

3-63 

1375 

3496 

■ 7 

Shuttle 

1 .81 

2308 

2433 

6 

Shuttle 

1 .81 

2308 

2433 

15 

Shuttle 

3-56 

1726 

2864 
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TABLE IV.- CHEMICAL SCHEME* 


Reaction 


1 . 

OH 

+ 

H 2 



= 

h 2 o 

+ 

H 

2. 

H 

+ 

OH 

+ 

M 

= 

HpO 

+ 

M 

3. 

H 

+ 

HC1 



r 

Cl 

+ 

H 2 

4. 

CO 

+ 

HOp 




co 2 

+ 

OH 

5. 

H 

+ 

ci 2 

+ 

M 

— 

HCI 

+ 

M 

6. 

0 

+ 

Cl 

+ 

M 

= 

CIO 

+ 

M 

7. 

0 

+ 

CIO 



= 

Cl 

+ 

°2 

8. 

HC1 

+ 

OH 



= 

HpO 

+ 

Cl 

9. 

N 

+ 

NO 



= 

n 2 

+ 

0 

10. 

NOo 

+ 

H 



r 

NO 

+ 

OH 

11 . 

CIO 

+ 

H 



= 

HCI 

+ 

0 

12. 

Cl 

+ 

Cl 

+ 

M 

□ 

Clp 

+ 

M 

13- 

H 

+ 

HOp 




H ? 

+ 

°2 

14. 

NO 

+ 

°3 



= 

no 2 

+ 

°2 

15. 

CO 

+ 

orf 



= 

COp 

+ 

H 

16. 

O3 

+ 

0 



= 

°2 

+ 

°2 

17. 

Cl 

+ 

°3 



= 

CIO 

+ 

°2 

18. 

0 

+ 

HCI 



r 

Cl 

+ 

OH 

19. 

H 

+ 

ci 2 



□ 

HCI 

+ 

Cl 

20. 

H 2 

+ 

HOp 



= 

HpO 

+ 

OH 

21 . 

NO 

+ 

CIO 



□ 

Cl 

+ 

NO. 

22. 

OH 

+ 

H0 2 



= 

°2 

+ 

Hp( 

23. 

O 


HO 2 



= 

OH 

+ 

°2 

24. 

NO 

+ 

O 

+ 

M 

- 

NOp 

+ 

M 

25. 

0 

+ 

0 

+ 

M 

= 

°2 

+ 

M 

26. 

OH 

+ 

OH 



— 

HpO 

+ 

0 

27. 

0 

+ 

H 2 



r 

OH 

+ 

H 

28. 

H 

+ 

O2 



= 

OH 

+ 

0 

29. 

0 

+ 

H 

+ 

M 

r 

OH 

+ 

M 

30. 

CO 

+ 

0 

+ 

M 

= 

COp 

+ 

M 

31 . 

H 

+ 

H 

+ 

M 

= 

h 2 

+ 

M 

32. 

CIO 

+ 

OH 



= 

ho 2 

+ 

Cl 

33- 

H 

+ 

°2 

+ 

M 

= 

HO 2 

+ 

M 

34. 

Cl 

+ 

H0 2 




HCI 

+ 

°2 

35. 

H 

+ 

HOp 



= 

OH 

+ 

OH 

36. 

N 

+ 

Op 




NO 

+ 

0 


Rate 


4 x IO -11 exp(-5500/RT) 

(1 x io _28 )/T 

8.8 x io -11 exp(-4622/RT) 

1 x io -11 exp(-3968/RT) 

(3 x 10 -29 )/T 
(1 x 10 _3 °)/T 
5.25 x io -11 exp(-992/RT) 

7.2 x io -12 exp(-3250/RT) 

2.7 x io - ]'! 

5.8 x io -10 exp(-l468/RT) 

3.0 x io -11 

1 .6 x IO -33 exp( 1600/RT) 

1 x io -10 exp( -1984/RT ) 

1.5 x 10 -12 exp(-2637/RT) 

5 x IO -13 exp(-600/RT) 

2 x io -11 exp(-4600/RT) 

3 x io -11 

2.1 x IO -11 exp( -7100/RT) 

2.09 x 10“ 10 

1 .0 x io -11 exp(-23 808/RT) 

1.7 x IQ" 11 

2 x io" 11 

3 x io -11 exp(-4960/RT) 

4.17 x IO’ 33 exp( 1860/RT) 

(1 x io~ 29 )/T 

1 x io -11 exp(-1000/RT) 

3 x IO -11 exp( -8 200/RT ) 

3 x io“ 10 exp(-16 500/RT ) 

(1 x io _29 )/t 

(2 x 10" 29 )/T exp(-4000/RT) 

(1 x 10 _29 )/T 

4 x io -11 exp(-2777/RT) 

4 x 10~ 32 exp( 992/RT ) 

1 x 10“ 10 exp(-1984/RT) 

3 x io~ 10 exp(-1984/RT) 

(1.1 x 10 _1 ^)/T exp(-6250/RT) 


A1C1, AlpO^, and AICI2 are included in the plume as inert species. Reac- 
tion rate coefficients for reactions 1, 2, 5, 15, 26, 27, 28, 29, 30, and 31 are 
from reference 8; reactions 3, 8, 10, and 25 are from reference 13; reactions 4, 
7, 11, 13, 20, 22, 23, 32, 33, 34, and 35 are from reference 9; reactions 9, 14, 
24, and 36 are from reference 11; reaction 16 is from reference 10; reaction 17 
is from reference 12; reactions 18, 19, and 21 were furnished by H. Hoshizaki of 
Lockheed Aircraft Corporation. Rates for reactions 6 and 12 are within the error 
bounds set by reference 9 but are lower than the forward rate recommended. 
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TABLE V.- MASS FLOW OF Al^ 


[The following table gives an indication of how closely the law of mass 
conservation is obeyed . AlgO^ is a chemically inert species in this 
model . The mass flow at the exit plane should equal the mass flow 
1 km downstream . ] 


Shuttle : 

Sea level (0.7 km) 

6-km altitude 
10-km altitude 
15-km altitude 

Cl enhanced - 10-km altitude 
Cl suppressed - 10-km altitude 
0 enhanced - 10-km altitude 
0 suppressed - 10-km altitude 
Combustion enhanced - 10-km altitude 
Ting-Libby model - 10-km altitude 

Titan: 

Sea level (0.7 km) 

18-km altitude 


Mass flow at 

Mass flow at 

exit plane , 

1 km downstream, 

g/sec 

g/sec 

9.938 x 105 

9.699 * 10 5 

9.890 x 105 

9.209 x 105 

1.141 x 10 ° 

1.164 x 10 ° 

1.240 x 10 6 

1.245 x 10*? 

1.141 x 10*? 

1.089 x 10^ 

1.141 * 10 6 

1.086 x 10° 

1.141 x 10*? 

1.026 x 10 ” 

1.141 * 10 6 

1.053 x 10° 

O V 
O 

X 

-=r 

1.118 x 10° 

1.141 x 10° 

1.015 x 10° 

5.204 x 10 5 

4.989 x 10 5 

4.905 x 10 5 

4.635 x 10 5 


TABLE VI.- REACTION RATE CHANGES INVOLVING ATOMIC CHLORINE 
[Reaction rate coefficient = CT n exp(B/RT)] 


Reaction 


Cl accelerated C 


1 . 

0 

+ 

HC1 



= 

Cl 

+ 

OH 

2.1 

X 

10- 11 

1 .7 

X 

10 -11 

2. 

Cl 

+ 

°3 



= 

CIO 

+ 

°2 

1.59 

X 

10-11 

3.09 

X 

10 -11 

3- 

HC1 

+ 

OH 



— 

H o 0 

+ 

Cl 

1 .0 

X 

10 '] 

4.0 

X 

10- ? 

4. 

H 

+ 

HC1 



— 

& 

+ 

h 2 

8.8 

X 

10-1° 

8.8 

X 

10 -12 

5. 

H 

+ 

Cl 

+ 

M 

= 

HC1 

+ 

M 2 

1 .0 

X 

10 "I,] 

1.0 

X 

loil 

6. 

Cl 

+ 

Cl 

+ 

M 

= 

Ci- 

+ 

M 

5.0 

X 

10 -34 

1 

X 

10" 32 

7. 

NO 

+ 

CIO 



— 

2 

Cl 

+ 

N0 ? 

1 .9 

X 

10- 11 

1 .5 

X 

10- 11 

8. 

Cl 

+ 

H0„ 



= 

HC1 

+ 

°2 

2.0 

X 

1°"J 2 

1 .0 

X 

10 "?Q 

9. 

0 

+ 

ci 2 

+ 

M 

r 

CIO 

+ 

M 2 

6.6 

X 

10 -31 

3-0 

X 

10 -29 

10. 

H 

+ 

Cl~ 



= 

HC1 

+ 

Cl 

2.09 

X 

i°"!° 

1 .28 

X 

10-JO 

11 . 

CIO 


OH 2 



— 

H0 5 

+ 

Cl 

4.0 

X 

10" ° 

1.33 

X 

10 -12 

12. 

0 

+ 

CIO 



= 

ci 2 

+ 

°2 

5.0 

X 

10 -10 

5.0 

X 

10 -12 


Cl slowed C 
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TABLE VII. 


REACTION RATE CHANGES 

INVOLVING ATOMIC OXYGEN 





Reaction 




C for 0 

accelerated 

C for C 

slowed 

1 . 

°3 

+ 

0 

= 

^2 

+ 

°2 

1 .0 

X 

10-!! 

4.0 

X 

10" 11 

2. 

NO 

+ 

0 + M 

“ 

NO, 

+ 

M 2 

2.0 

X 

10-33 

8.0 

X 

10-33 

3. 

0 

+ 

0 + M 

= 

V 

+ 

M 

3-0 

X 

10" 31 

3-0 

X 

10 -28 

4. 

OH 

+ 

OH 

= 

H,0 

+ 

0 

5.0 

X 

10" 11 

2.0 

X 

10“ 12 

5. 

0 

+ 

h 2 

= 

OH 

+ 

H 

9.7 

X 

10" 12 

8.7 

X 

10-11 

6. 

H 

+ 

°2 

= 

OH 

+ 

0 

1.1 

X 

io i n 

1 .2 

X 

10" 1 ? 

7. 

0 

+ 

H + M 

r 

OH 

+ 

M 

3-0 

X 

10 -29 

3.0 

X 

10" 28 

8. 

N 

+ 

NO 


n 2 

+ 

0 

5.07 

X 

10 ~]l 

1.5 

X 

10"11 

9. 

N 

+ 

°2 

- 

NO 

+ 

0 

1 .4 

X 

10" 14 

8.3 

X 

LO 

— 

1 

O 

10. 

0 

+ 

HC1 

= 

Cl 

+ 

OH 

1 .7 

X 

10- 11 

2.2 

X 

10-11 

11 . 

CIO 

+ 

H 

= 

HC1 

+ 

0 

3-0 

X 

10" 10 

3.0 

X 

10-13 

12. 

0 

+ 

CIO 

— 

Cl 

+ 

°2 

5.25 

X 

10" 12 

5.0 

X 

io-i° 

13- 

0 

+ 

H°2 

□ 

OH 

+ 

°2 

1 .0 

X 

10" 12 

9.0 

X 

10-10 

14. 

CO 

+ 

0 + M 

— 

co 2 

+ 

M 2 

3.0 

X 

1 0“3 1 

8.0 

X 

C\l 

O 

15. 

0 

+ 

Cl + M 

— 

CIO 

+ 

M 

6.66 

X 

10-31 

3.0 

X 

10 -29 


TABLE VIII.- REACTION RATE CHANGES INVOLVING COMBUSTION 

[All other reactions are the same as in table IV. 
Rate = CT n exp(B/RT)] 





Reaction 





C for 

combustion 











accelerated 


1 . 

OH 

+ 

H 2 


- 

HpO 

+ 

H 

4.0 

X 

10 '!! 


2. 

0H 

+ 

OH 


= 

H,0 


0 

5.0 

X 

10*11 


3. 

0 

+ 

H 2 


- 

oi 

+ 

H 

8.7 

X 

IQ -! 1 


4. 

H 

+ 

°2 


= 

OH 

+ 

0 

1 . 1 

X 

10" 9 


5- 

H 

+ 

H 

+ M 

— 

H, 

+ 

M 

1 .5 

X 

10-2 


6. 

H 

+ 

O 

+ M 

= 

oi 

+ 

M 

3-0 

X 

1 0 -28 


7. 

H 

+ 

OH 

+ M 

= 

HpO 

+ 

M 

2.0 

X 

io iJ 


8. 

0 

+ 

0 

+ M 

= 

°2 

+ 

M 

3-0 

X 

10" 28 


9. 

CO 

+ 

OM 


- 

COp 

+ 

M 

5-4 

X 

10- 


10. 

CO 

+ 

0 

+ M 

= 

COp 

+ 

M 

8.0 

X 

10“ 29 


. 
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TABLE IX.- N0 X PRODUCTION (AT 1 km DOWNSTREAM FROM EXIT) 



NO, g/m 

N0 2 , g/m 

Shuttle : 



Sea level (0.7 km) 

1052 

12.88 

6-km altitude 

116 

1.5 

10-km altitude 

25 

.48 

15-km altitude 

6 

.14 

Cl enhanced - 10-km altitude 

29 

.166 

Cl suppressed - 10-km altitude 

37 

.82 

0 enhanced - 10-km altitude 

30 

.82 

0 suppressed - 10-km altitude 

20 


Combustion enhanced - 10-km altitude 

32 

.64 

Ting-Libby model - 10-km altitude 

41 

.20 

Titan: 



Sea level (0.7 km) 

489.6 

6.34 

18-km altitude 

21 .7 

. 12 


TABLE X.- CHLORINE PARTITIONING PERCENT BY WEIGHT OF REACTING COMPOUNDS 
CONTAINING CHLORIDE AT 1 km DOWNSTREAM 



HC1 

C1 2 

Cl 

CIO 

Shuttle : 





Sea level (0.7 km) 

89.75 

10.10 

0.14 

0.01 

6-km altitude 

88.69 

10.68 

.62 

.01 

10-km altitude 

85.54 

12.65 

1 .80 

.01 

15-km altitude 

78.06 

15.13 

6 .80 

.11 

Cl enhanced* - 10-km altitude 

83.91 

11.71 

4.37 

.02 

Cl suppressed* - 10-km altitude 

84.67 

14.70 

.62 

.03 

0 enhanced - 10-km altitude 

86.74 

11.22 

2.02 

.01 

0 suppressed* - 10-km altitude 

84.34 

13-93 

2.14 

0 

Combustion enhanced* - 10-km altitude 

84.55 

13-27 

2.17 

.01 

Ting-Libby model - 10-km altitude 

92.49 

7.11 

.45 

0 

Titan : 





Sea level (0.7 km) 

89.90 

9-90 

.20 

0 

19-km altitude 

74.56 

12.35 

13-12 

.01 


* 

Plume reactions are still occurring, these numbers will change further 


downstream . 
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TABLE XI.- CARBON PARTITIONING 


Shuttle : 

Sea level (0.7 km) 

6-km altitude 
10-km altitude 
15-km altitude 

Cl enhanced - 10-km altitude 
Cl suppressed - 10-km altitude 
0 enhanced - 10-km altitude 
0 suppressed - 10-km altitude 
Combustion enhanced - 10-km altitude 
Ting-Libby model - 10-km altitude 

Titan : 

Sea level (0.7 km) 

18-km altitude 


CO at 

C0 2 at 

CO 


exit 

exit 

(x = 1 km) , 

CO 2 > 

plane , 
g/m 

plane , 
g/m 

g/m 

g/m 

356 

51 

50.5 

29 828 

334 

48 

17 

4 371 

312 

45 

39-5 

3 283 

319 

45 

37 

2 705 

312 

45 

22 

3 285 

312 

45 

30 

3 262 

312 

45 

18 

3 083 

312 

45 

21 

3 176 

312 

45 

26 

3 358 

312 

45 

13 

2 942 

213 

25 

38 

17 698 

173 

21 

18 

1 212 
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TABLE XII.- AMOUNT OF COMPOUNDS 


Amount of certain compounds 1 km downstream of exit plane, g/m 


Titan , 

Titan , 

Shuttle , 

Shuttle , 

Shuttle , 

Shuttle , 

0.9 km 

18 km 

0.7 km 

6 km 

10 km 

15 km 


HC1 

7 565 

Cl 2 

833 

CIO 

.0 

Cl 

16.7 

NO 

489 

no 9 

6.3 

TO 

38.3 

co 2 

‘ 17 698 

u 1 2°3 

11 194 

h 2 

.0 

h 2 o 

11 791 

0 

.0 

OH 

.0 

j 


(Still reacting) 


708 

1' 

543 


.5 


21 .6 


954 

l 

12.9 


50 


828 

4 

807 

3. 


(Still reacting) 















TABLE XII.- Concluded 



Amount of certain compounds 1 km downstream of the exit plane, g/m 
for 10-km shuttle with parameter variations 


Cl reactions 
faster 

Cl reactions 
slower 

0 reactions 
faster 

0 reactions 
slower 

Combustion 

faster 

Ting-Libby 

HC1 

1419 

1428 

1376 

1379 

1465 

1391 

£s 

198 

248 

178 

221 

230 

107 

.04 

.05 

.15 

.04 

.18 

.002 

Cl 

74 

10.5 

32 

35 

37.6 

6.7 

NO 

29 

37 

30 

20 

32 

40.9 

N0„ 

.16 

.8 

.8 

.17 

.64 

• 19 

co 2 

22 

30 

17 

21 

26.5 

12.8 

CO 

3285 

3262 

3084 

3176 

3358 

2942 

Al 2°3 

h 2 

2420 

2412 

2268 

2339 

2478 

2146 

.002 

.008 

.000 

.000 

.005 

.000 

H,0 

2324 

2314 

2167 

2244 

2377 

2029 

0 

.01 

.001 

.005 

.002 

.007 

.000 

OH 

• 03 

.003 

.009 

.010 

.014 

.000 

(Still reacting) 





Mass Flow of 


CO 


o 

<M 


<1 



( 


Figure 2.- Predicted mass flow 
distance in plume of a shut 
unmodified and modified LAP 






Quantity of A1 2 C>3, g/i 



Percent deviation from center line values 




4 

bJD 

of 

<x> 

•pH 

o 

& 

CO 


O 



0 40 80 120 160 200 240 280 
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Figure 4.- Production of N0 X as a function of downstream distance in plume 
of a shuttle motor at an altitude of 10 km by using Donaldson-Gray and 
and Ting-Libby eddy-viscosity models. 





Mole fraction of species along center line 



Figure 8.- Mole fractions along center line of several species related to 
combustion processes as a function of downstream distance in plume of 
a Titan motor at 0.9 km. 
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Quantity of species, 






Figure 11 Quantity of several species related to combustion processes as a 
function of downstream distance in plume of a shuttle motor at altitudes 
of 0.7 km and 15 km. 




0 40 80 120 160 200 240 280 

x/r 


Figure 12.- Quantity of several species containing chlorine as a function of 
downstream distance in plume of a shuttle motor at altitudes of 0.7 km 
and 15 km. 





x/r 

Figure 14.- Mole fractions along center line of several species related to 
combustion processes as a function of downstream distance in plume of a 
shuttle motor at altitudes of 0.7 km and 15 km. 
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Figure 15.- Mole fractions along center line of several species containing 
chlorine as a function of downstream distance in plume of a shuttle 
motor at altitudes of 0.7 km and 15 km. 




x/r 

Figure 16.- Mole fractions along center line of several species related to NO 
as a function of downstream distance in plume of a shuttle motor at alti- 
tudes of 0.7 km and 15 km. 



Velocity, m/sec 


3 2 1 0 1 2 3 45 30 15 0 15 30 45 

y/r y/r 



Figure 17.- Velocity profiles at several downstream stations in plume of a 
shuttle motor at an altitude of 0.7 km. 
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Figure 18.- Velocity profiles at several downstream stations in plume of a 

shuttle motor at an altitude of 15 km. 



Temperature, K Temperature, 


2500 



Figure 19 •- Temperature profiles at several downstream stations in plume of 
a shuttle motor at an altitude of 0.7 km. 
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